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Abstract 

Following an idea of Rowland [Row] we give a conjectural way to gen- 
erate increasing sequences of primes using "gcd-algorithms". These al- 
gorithms seem not so useless for searching primes since it appears we 
found sometime primes much more greater than the number of required 
iterations. In an other hand we propose new formulations of famous con- 
jectures from the additive theory of numbers (the weak twin prime conjec- 
ture, the Polignac conjecture, the Goldbach conjecture or the very general 
Schinzel's hypothesis H). For the moment these are experimental results 
obtained using pari-gp |PAR| . 

Introduction 

In |Row| the author solved a question related to the sequence defined recursively 
by i?(l) = 7 and forn > 2 by: 

• R{n) = R{n - 1) + gcd(i?(n - 1), n) 

Namely he showed that R{n) — R{n — 1) is always a prime number or 1. At first 
glance this recursion is not very useful for finding primes since these differences 
generate primes quite randomly and the prime values are less than the number 
of required iterations. However, looking closer to Rowland recursion we found a 
way to exhibit the records values in the sequence of differences R{n) — R{n — 1) 
(these differences are given by yll37613 in [Slo ). Namely n + 1 when R{n) = 
2n + 2 are exaclty those records values. So this gives a method to exhibit 
primes in increasing order (and the growth is exponential). Those n + 1 such 
that R{n) = 2n + 2 begin: 

5, 11, 23, 47, 101, 233, 467, 941, 1889, 3779, 7559, 15131, 30323, 60647, . . . 

Although this seems provable using Rowland reasoning, it looks like it is un- 
known since it is not mentioned before in |Row| and this sequence of records 
values is not in |Slo| . This is clearly a nice fact and we try to generalise this 
observation. In order to facilitate the computations we use the absolute value 
and run the algorithm "backwards" as we will see. So we succeeded to extend 
somewhat this result with the conjectures 1 and 2 and propose a method for 
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building primes (section 1). In a conjecture 3 (section 2) we claim there are 
sequences of lesser of twin primes growing very fast. Then we provide a con- 
jecture 3bis (cf. 2.4) related to the Polignac conjecture and a conjecture Ster 
(cf. 2.5) dealing with the prime triplet conjecture. In a conjecture 4 (section 
3) we propose a way to find very big primes compared to the number of itera- 
tions which could improve the conjecture 2 and be an efficient tool for searching 
primes. In a conjecture 5 (section 4) we propose a first constructive way to 
prove the Goldbach conjecture but with some restiction due to an exceptional 
set of non working cases. In a conjecture 6 (section 5) we propose to summarize 
all our observations with a reformulation of the Schinzel's hypothesis H. The 
conjecture 7 (cf. 5.2) gives another approach of the conjecture 6. 

We also discuss the Shevelev conjecture who extends Rowland idea for twin 
primes. Indeed, V. Shevelev |She| introduced the related sequence (^166944 in 

m) ■■ 

• S{1) = 2 

• S{n) = S{n - 1) + gcd(S'(7i - 1), n - 1 + (-1)") 

And noticed that the records values of differences (^166945 in |Slo| ) : 
7,13,43,139,313,661,1321,... 

gives for terms > 7 the greater prime of some twin primes pairs . As for Rowland 
sequence it can be seen that the indices where these records occur is given by the 
sequence of 1 such that S{n) = 2n + l (except for 7) and (n, S{n) — S{n— 1)) 
is then a twin pair of primes. Using this algorithm backwards the conjecture 8 
(section 6) propose a sligthly different way than the conjecture 3 or 6 or Shevelev 
conjectures to prove the weak twin prime conjecture. 

Although Shevelev managed to do it, it is not easy to generalise this observa- 
tion using the "forward" original Rowland recursion since we need to find where 
the records occur. However using the absolute value we shall see it is easy to 
obtain many increasing sequences of twin primes or of triplet of primes since we 
just have to check indices where zeros occur. Thus in a conjecture 9 (section 7) 
we propose a method for generating primes m-uplet of any type using periodic 
sequences in the gcd algorithm. Then we extend the idea to a family of poly- 
nomials giving another version of the Schinzel's hypothesis H. This allows us to 
perform easier computations than using the conjecture 6 for searching m-uplet. 

We finally merge our "backward" recursion and Shevelev idea for approaching 
the Goldbach conjecture in the conjecture 10 (section 8) in a nicer way than the 
conjecture 5. There is apparently no more exception for n large enough. This 
right way to deal with this old conjecture using gcd-recursion was very hard to 
find. We also discuss a less known but hard conjecture of Legendre (section 9). 

There is a curious fact in this study. It appears our methods work often 
better for large integer since once we found a good starting value we generate 
only primes, despite the probability to pick up big primes among large integers 
goes to zero. 
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1 Variation on Rowland recursion 



After some attempts we arrive to this recursion seeming generating primes in a 
general way starting always with the same initial value 1. For a given integer 
value m > 1 we define the sequence (a(n))^>j^ by a(l) = 1 and for n > 2 by the 
recursion: 

• a{n) = \a{n — 1) — gcd {a{n — 1), mn — 1)|. 
Now let us consider the values of n such that we get: 

• a(n) = ^ a(n + 1) = mn + m — 1. 

We claim that Vm > 1 this sequence of indices n gives rise to an infinite sequence 

(^'m(fc))fe>i. 

Example 

For TO = 1 the scquence(a(n))„>j begins: 

1, 0, 2, 1, 0, 5, 4, 3, 2, 1, 0, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1, 0, 23, 22, 21, 20, 19, . . . 
And the sequence (6i {k))^.^^ of indices where zeros appear in (o(n))„>i begins: 
2, 5, 11, 23, 47, 79, 157, 313, 619, 1237, 2473, 4909, 9817, 19603, 39199, 78193, . . . 
These listed numbers are prime numbers. This led us to a first conjecture. 

1.1 Conjecture 1 

We claim: 

• bi{k) is prime for A; > 1 and bi{k) ~ c2'' {k oo) with c = 1.186 . . . 

See the APPENDIX 1 for a table supporting this conjecture. Although this 
conjecture is similar to the record sequence mentioned in the introduction, it is 
more interesting to us. Indeed we were able to generalise the result and much 
more seems true. So we make a stronger conjecture. 

1.2 Conjecture 2 

In general wc claim that for any to > 1: 

• mbm (A;) + TO — 1 is prime for k large enough (usually fc > 2 is working for 
small m). 

• mbm{k) + m — 1 Cmim + 1)'' (fc — >■ oo) with Cm > 0. 



3 



Additional examples supporting the conjecture 2 

For m = 3 the sequence (a(Ti))„>i begins: 

1, 0, 8, 7, 0, 17, 16, 15, 14, 13, 12, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1, 0, 71, . . . 
And the sequence {b3{k)) ^.-^-^ begins: 
2, 5, 23, 89, 337, 1335, 5307, . . . 

Next the values 363 (k) + 2 appear to be prime values for fc > 2: 
8, 17, 71, 269, 1013, 4007, 15923, 63521, 253949, 1014317, . . . 
And 363(71) + 2 ~ C34" (n 00) with C3 = 0.96 . . . 
For m = 4 the sequence (a(«))„>i begins: 

I, 0, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1, 0, 59, 58, 57, . . . 
And the sequence {b4{k))^y-^ begins: 

2, 14, 62, 314, 1574, 7846, 38020, . . . 

Next the values 464 (/c) + 3 appear to be prime values only for k > 1: 

II, 59, 251, 1259, 6299, 31387, 152083, 758971, 3790651, 18953251, . . . 
And 464(71) + 3 ~ C45" (n -> 00) with a = 1.9408 . . . 

1.3 An efficient algorithm for finding primes? 

Note we find sometime primes greater than n after making n iterations. So 

it could be an efficient method for finding primes since the gcd algorithm is 
well known and "fast" computation can be performed. For instance if w = 28 
we compute 1000000 terms of the sequence (a(T^))„>i • In this range a(n) 
vanishes 5 times. This allows us to compute (628('^))i<„<5 wich gives 5 values 

of 28628 (n) + 27: 

83, 1147, 31891, 924811, 26819491. 

All these values are primes and the last one 26819491 gives a prime number 
larger than the 1000000 iterations (see APPENDIX 2 for more exemples). More- 
over it appears this kind of algorithm can be adapted for finding bigger primes 
as shown thereafter. 
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1.3.1 A simple rule of construction 

Observe from our definition we have for n > 2 (letting b = bm )'■ 

. a{b{n)+l) = (m+l)a(6(n-l)+l)+m+mE?L"i;'_i)+i {a{j + 1) - a{j) + 1) 

Hence we have an "almost" recurrence relation between two consecutive records 
values which are conjectured to be prime for n large enough. This formula 
explains also why these record values are growing like (m + 1)" since it appears 
we usually don't need to compute all terms in the sum. Indeed experiments show 
that a(j + 1) — a{j) = —1 for b{n — 1) < tiq < j < b{n) — 1 and this ng depends 
on b{n — 1) and stays "near" from this value, i.e. we claim that no — b{n — 1) <C 
\/b{n — 1) and is often much more smaller. So we can launch a computation 
and when the computer returns a{j + 1) — a{i) = —1 sufficiently "often" we may 
stop the computation and suspect we are beyond this no- Therefore we built 
perhaps a prime greater than a starting prime. It could then be interesting 
to know rules in order to choose "good" values of m making no the closest of 
6(n — 1) as possible. 

1.3.2 Exemples 

For instance consider rn = 10 and the sequence of a„ — a„_i when |a„ — a„_i| > 
1 wich yields somewhere for a„ — a„_i: 

• 43213789, -3, -13, -15241, -43, -1889, -3, -433, -113, -3, -5827, -247 

The positive terms are our record values a(b{k) + 1) and are primes values. 
Here the computation returns nothing more using 10^ iterations. Hence we are 
perhaps beyond our no and the next record would be given by 11 x 43213789 + 
10-10x(2 + 12 + 15240 + 42 + 1888 + 2 + 432 + 1 12 + 2 + 5826 + 246) = 4751 13649 
which is a prime number and, as expected, our next record value. So we have 
found a prime number 11 times greater than the given prime 43213789 with few 
iterations. 

But let us see this with a more striking exemple and take ni = 100000. We 
compute |a„ — a„_i| > 1 for n < 10^ wich yields for these a„ — a„_i: 

• 299999, -59, 29994499999, -3, -7, -53, -3 

So we suspect we have nothing more after —3 until the next record value. Thus 
100001 X 29994499999 + 100000 - 100000 x (2 + 6 + 52 + 2) = 2999479988299999 
should be the next record and it is indeed a prime value. So with few effort 
we found a prime number 100001 times greater than another one. A simple 
routine under pari-gp found easily big primes with this method (see the end of 
the APPENDIX 2). 

1.3.3 Probability to got a prime value from the first record 

We also observe something which could have practical use for searching primes 
if people work together like for the GIMPS where they share computer power. 
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We keep previous notations and our definition of a„ (depending on m) so that 
the first record is simply 3m — 1. Then we tried to estimate the chance to got 
a prime value with the second record value for various m after making [m"J 
iterations and where < a < 1 . i.e. for a given m we define: 

• fa ifn) = 1 if the second record value which equals (m + 1) (3to — 1) + to + 
TO^j™^^ (a{j + 1) — a{j) + 1) is a prime value and otherwise ,fa{ni) = 0. 

• L{a) = lim„^oo ^ ELi /"(^) 
Then we conjecture that: 

• L{a) exists with simply L{a) =^ a and so we know the probability to got 
a prime with this method. 

To see this last point with experiments suppose E \& a, set of 30000 random 
values TO satisfying lO^'' < to < 10^^ . We take a = ^ so that we have to 
make around 100 iterations only to got the value (to + 1)(3to — 1) + to + 
'Ti^l™^^ {a{j + 1) — a{j) + 1). We then launch the computation 3 times (so 
that E changes each time) and we plot the different graph of ^ X]r=i fa{'nii) 
(graphs are black, blue and green) where toi < TO2 < ... G E compared to the 
graph of y = i (red). 

(fig-1) 



0,15 




o.i: 



□ ,13 



This supports the claim L{a) = a. In other word this means that if we take 
[a~^J random values of m we are fairly certain the formula (to + 1)(3to — 1) + 

TO + TO^I™^^ (a(j + 1) — a(j) + 1) will produce at least a prime number among 
them. Note we may start form the second record value (or the third...) instead 
of the first but we aim to work with big values of to so it seems inappropriate 
to compute the second record. 
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Remcirk 

Although this method for finding primes seems satisfying it is not evident to see 
whether it is very efficient since we need to combine two computations (one for 
finding the candidates and one for testing primalty) and so far we find no rule 
in order to force uq to stay very close from b{n — 1). Perhaps the conjecture 4 
would be better for that purpose (see 3.4.). By the way this kind of algorithms is 
interesting in its own since variations of the Rowland-Shevelev algorithm led us 
to formulate in a new way old conjectures from the additive theory of numbers. 
Hereafter we come across conjectures like the Polignac conjecture, the Schintzel 
hypothesis H or the Goldbach conjecture and our approach suggests analytic or 
probabilistic study. 

2 An increasing sequence of twin primes 

This was somewhat surprising to find such a sequence since it is not known 
whether there are infinitely many twin primes. Here we consider the recursion 
a(l) = 1 and for n > 2: 

• a(n) = \a{n — 1) — gcd(a(n — 1), — 1)| 

2.1 Conjecture 3 

Suppose p > 2 is a prime number then we claim: 

• there are infinitely many values of n such that a{n) = 0. 

• for n large enough (usually n > 2 is working for small p) we have a{n) = 
n is prime and — is prime. 

2.2 Corollary 

There are infinitely many twin primes since the conjecture yields for p = 2 : 

• for n large enough we have a{n) = n is prime and n + 2 is prime 
and this happens infinitely many times. 

2.3 Tables (in the sequel 5{n) = 1 if n is prime and otherwise.) 

p = 2 



Values of n such that a{n) = 


n + 2 


S{n) 


6{n + 2) 


2 


4 


1 





11 


13 


1 


1 


137 


139 


1 


1 


19181 


19183 


1 


1 


367953497 


367953499 


1 


1 
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p = 3 



Values of n such that a(n) = 


n 


S{n) 




2 


23 


1 


1 


23 


601 


1 


1 


119333 


142432291 


1 


1 



We can also start from another initial value. For instance let p = 2 and 
choose a(l) = 2 this gives 



Values of n such that a{n) = 


n + 2 


S{n) 


S{n + 2) 


3 


5 


1 


1 


17 


19 


1 


1 


281 


283 


1 


1 


79559 


79561 


1 


1 


6329815697 


6329815699 


1 


1 



See the APPENDIX 3 for more experiments supporting the conjecture for 
p = 2 and various starting values. 

2.4 Conjecture 3bis 

In the same vein we find an algorithm seeming generating primes pairs of type 
{p,p + 2m). Let a(l) = Am? and define: 

• a(n) = \a{n — 1) — gcd(a(n — 1), n(n + 2m))| 
Then we claim that for n large enough: 

• a(n) = + 1 and n + 2m + 1 are primes and this happens infinitely 
many times and there are infinitely many pairs (n + l,n + 2m + 1) of 
consecutive primes. 

Thus the Polignac conjecture would be true. Here a table for 2m = 4 



Values of /) wIkui (i\n) = 




()(" + 1) 


d(ii + 5) 


12 


17 


1 


1 


192 


197 


1 


1 


38196 


38201 


1 


1 


1459118862 


1459118867 


1 


1 



See APPENDIX 4 for experiments with other small values of m. The con- 
jecture 9 will give an easier way for computation of such pairs of primes since 
we avoid the square in the gcd. 
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Remark It is also possible to use Rowland-Shevelev recursion for generating 
increasing sequences of twin primes or things like that. For instance let; 

• ai = 2 and a„ = a„_i + gcd(a„_i, n(n — 2)) 

Then the sequence of records for the differences a„ — a„_i yields an increasing 
sequence of lower of twin primes. See the end of the APPENDIX 3. 

2.5 Conjecture 3ter 

We generate primes triplet of type {p,p + 2,p + 6). Let a(l) = 4 and define: 

• a(n) — \a{n — 1) — gcd(a(n — 1), n(n + 2)(n + 6))| 
Then we claim 

• a{n) = =>(n + 1, n + 3, n + 7) is a prime triplet. 
Here a table 



Values of n such that a{n) =0 


S{n+1) 


S{n + 3) 


S{n + 7) 


4 


1 


1 


1 


40 


1 


1 


1 


82006 


1 


1 


1 



Although it is hard to perform convincing experiments, it is coherent with 
our previous conjectures and some very general rule should exist. This conjec- 
ture can be extended to any sort of prime triplet and to r?>uplet. The conjecture 
9 is better fitted for practical computation of m-uplet. 

3 Finding big primes 

Here we merge the conjectures 2 and 3 to obtain another conjectural way to un- 
earth increasing sequences of primes. The rate of growth is multiple exponential 
and this could produce very big primes compared to the number of required it- 
erations. We tlien discuss probabilities issues and propose a method to got very 
big primes similarly as what is described in 1.3. 

3.1 Conjecture 4 

Suppose p > 2 is prime and let a(l) = p and: 

• a{n) = \a{n — 1) — gcd(a(n — l),pn^ — 1)| 
Then we claim: 

• there are infinitely many values of n such that a{n) = 0. 

• for n large enough (usually n > 2 is working for small p) we have a{n) = 

=>p{n + 1)^ — 1 is prime. 

See APPENDIX 5 for the begining of tables for some values of p. Could we by 
chance get so many prime numbers? 
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3.2 Why stopping here? 

In fact it seems one can go further. Suppose a(l) = 2 and: 

• a(n) = |a(n — 1) — gcd(a(n — 1), 2n^ — 1)| 
Then we claim: 

• there are infinitely many values of n such that a{n) = 0. 

• o(n) = =>2(n + 1)^ — 1 is prime for n large enough. 
Here the begining of the table 



Values of n such that a{n) =0 


2(n+ 1)^ - 1 


<5(2(n + 1)^-1) 


3 


127 


1 


125 


4000751 


1 


4000877 


128084306502569672303 


1 



However it is not easy to generalise this one. We can also find good initial 
values and/or a good factor before n^. For instance let w(l) = 3 and: 

• a(n) = |a(n - 1) - gcd(a(n - 1), lOn^ - 1) | 

Then we get 



Values of n such that a{n) = 


10(n+ 1)^ - 1 


(5(10(n+ 1)^ - 1) 


4 


1249 


1 


1240 


19112405209 


1 



And we can provide this other impressive example with exponent 7. Suppose 
a(l) = 3 and: 

• a(n) = |a(n — 1) — gcd(a(n — 1), 2n'^ — 1)| 



Values of n such that a{n) =0 


2(n + 1)^ - 1 


6{2{n+iy-l) 


2 


4373 


1 


4352 


59231218330987879606185473 


1 



3.3 How many chances we have to catch a prime? 

Continuing this way it seems possible to find very big primes using the {k, b, c) 
recursion: 

• a(l) = k and a{n) = \a{n — 1) — gcd(a(n — 1), bn'^ — 1)| 

with suitable choices of (fc, &, c). A natural question is then: what is the chance 
to get a prime when the algorithm reachs the first zero starting with any value 
A''? To evaluate this chance let rk = min {i > 1 | a{i) = 0} and: 
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. T{N) = i# {fc I 1 < A; < NkS{b{rk + lY - 1) = 1} 

wich represents the chance to get a prime reached by the algorithm starting with 
A^. As A'^ — >■ oo it appears this chance is not zero. For instance if (6, c) = (2, 2) 
we have T(7V) ~ 0.8 as iV — oo. Which gives an efficient method to got big 
primes since the first zero is reached after a number of iterations of order N and 
thus we have more than 80% of chance to got a prime of size 2N'^. Here is a 
graph supporting this claim. We plot T{N) for N = 1,2,3,..., 20000 and for 
(6,c)-(3,2) (pink) (6,c) = (5,2) (blue) 



(fig-2) 




And we believe limTv^oo T(iV) — L{c) exists and depends only on c and it 
is clear that c < c' L{c) > L{d). 

3.4 A rule of construction 

As seen before in 1.3. there is here also a general relation between consecutive 
records. Namely we still consider for any initial value: 

• a{n) = \a{n — 1) — gcd(a(n — 1), brf — 1)| 

and we define (w(/c))j,>-^ as the sequence of values taken by a„ — a„_i when 
|on — fln-il > 1- Then the records values of {w{k))j.-~,-^ are given by k such 
that w(A;) > 0. Let us now write this increasing sequence of k using a sequence 
('^j)j>i- Then we have the following simple relationship between 2 consecutive 
records values w{aj) and ^^{aj+i) (details ommited): 

• = b {w{a,) + 1 + (^^) + ^f^;-; (a,+i - a, + 1))" - 1 
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Thus as in 1.3.1. experiments show that it isn't necessary to compute all terms 
in the sum to got the next record value. Indeed there is again a value no 
conjectured to be closed of aj such that Uj < no < i < ctj+i Oj+i — Oi = — 1. 
Therefore this gives sometime an efficient method for building a bigger prime 
from a prime record value or a non prime record value. The quadratic case 
seems well working. For instance let us consider this quadratic case: 

• a(l) = 1 and a{n) = |a(n — 1) — gcd(a(n — 1), 32n^ — 1)| 

Then we get the sequence of values a„ — Un-i for those n < 10^ such that 

127, -7, -17, -7, -7, 294911, -1289, 2760686028799, -113, -103, -7, -7, -113. 

The 3 first records (127, 294911, 2760686028799)arc prime values and supposing 
there is no more value until the next record (say X) the formula above yields: 

X = 32 (^2760686028800 + ( 2760686028800 ^V2 _ ;^]^2 - 102 - 6 - 6 - 112^^ - 1 

Giving X = 243884447023448880167715967 which is still a prime value. We 
provide also an exemple for the cubic case: 

• a(l) = 2 and a{n) = \a{n — 1) — gcd(a(n — 1), 2n^ — 1)| 

Then we get the sequence of values a„ — a„_i for those n < 10^ such that 

-3, 53, -5, -3, 265301, -109, -31, -17, -3, -5, -3. 

So we suspect there is no more value until the next record value (say X) and 
the formula above yields: 

X = 2 (^265301 + 1 + ( 265301+1 ^ _ 108 - 16 - 2 - 4 - 2^ ^ - 1 

giving X = 37299785868725741 which indeed is the next record value and is a 

prime value. 

For higher exponent it seems less easy to find many working exemples but 
we think it would be worth to explore this method further in order to check 
its possible efficiency. The main question would be: are there any conditions 
forcing records values to be prime values and making no very close from the 
working record value? 

3.5 Generating big twin primes 

We can also do the same kind of task (cf. 3.3.) for twin primes borrowing from 

Shevelev the idea for a quadratic case. Let: 

• a(0) =k and a{n) = a{n - 1) - gcd(a(n - 1), 2n2 + (-1)") . 
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Define: 

• rt^,^n{N) - i# {fc | 1 < fc < NkS{2{rk + 1)^ - l)(5(2(rfe + l)^ + 1) = l} 

wich represents the chance to get a pair of twin primes reached by the algorithm 
starting with N. We plot below Ttwin{N) for = 1, 2, 3, ... , 20000 

(fig-3) 




And the graph in the range 2<N< 30000 of '^6{2{rN + 1)^ - l)(5(2(rjv 



(fig-4) 
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The fact the x axis is blue means there are many zero values. 

Roughly speaking it shows we have 50% chance to come across a twin prime 
of size N"^ after N iterations starting from a large random value of N . Certainly 
one should search for simple conditions on TV in order to increase this chance. 
For instance let us start from 2iV^ when 2iV^ — 1 and 27V^ + 1 are already primes. 
This sequence of N begins: 

• 3, 6, 21, 24, 36, 42, 45, 87, 102, 132, 153, 186, 204, 228, 237, 273, 297, 300, 321, ... 
Let now v{n) denotes the n-th term of this sequence and let: 

. T„(7V) - i# {fc I 1 < fc < iV&<5(2(r„(,) + 1)2 - l)(5(2(r,(fe) + if + 1) = l} 
We plot below T.„(iV) for = 1, 2, 3, ... , 140 



(fig-5) 




So it seems we have slightly more chances to got a twin prime pair starting 
with 2v{nf (~ 70% of chance to got a twin prime of order n^) than starting with 
n ('^ 50% of chance producing primes of order v?). Although some computation 
is needed to got the sequence v{n) it could be an efficient method for generating 
big primes. More importantly, this observation is a striking one regarding the 
conjectures like 2 or 3. Indeed, it appears that when we catch a prime with 
a given property and run the gcd-algortihm another time (starting around this 
value) we have more chance to get a new prime with this property. Heuristically 
this should explain why there are apparently infinite "chains" of primes generated 
by the algorithm in conjectures 2, 3 or similar ones when we have a good starting 
value. Repeating the process seems to force the algorithm to reach primes every 
time once we are on a right track. 
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4 Conjecture 5 : on the Goldbach conjecture 

We propose a first constructive way to prove this famous conjecture. This is 
somewhat unsatisfactory since there are very few non working starting values 
in the range where we performed the computation. Consider N > 2 and let 
a(l) = iV - 2 and define for n > 2: 

• a(n) = a{n - 1) - gcd(a(n - 1), {n - 1){2N -n + 1)) 

Then we claim there is always a unique gj^ g{2,3,...,A^ — 2} such that : 

• a(ffAr) = 
and we have: 

• (7 AT and 2N — are simultanuously primes except for very few iV ( a set 
conjectured to be of measure zero). 

Hereafter we plot S^5{gN)5{2N - qn) for 2 < iV < 30000. 

(fig-6) 
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We can see on the x axis some zero values which become very sparse when 
N increases (perhaps there is no more zero value for N sufficiently large). In 
the conjecture 10 we provide a variation of this conjecture where clearly there 
is no exceptional set of non working values for N large enough. 

5 The Schinzel's hypothesis H 

Finally this kind of method should have considerable application. One of them, 
regarding what we discuss before, could be a new version of the Schinzel's hy- 
pothesis H |Sch| . Indeed, suppose P{x) is a polynomial with integer coefficients 
and define the sequence S as follows. 
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• S{1) e N 



• S{n) = \S{n - 1) - gcd{S{n - 1), P(n))| 

Under some assumptions the values of n such that S{n) = imply P{n + 1) has 
nice arithmetical properties. This is what we have seen previously. But let us 
consider three additional examples. 

P{x) + 1 

We claim that for a suitable starting value S{1) such as 5(1) = 2 we have: 

• S{n) = for infinitely many values of n. 

S{n) = P{n + 1) is a prime number. 
Here a table with the first values 



Values of n such that S{n) = 


P{n + 1) 


S{P{n+l)) 


3 


17 


1 


13 


197 


1 


203 


41617 


1 


41813 


1748410597 


1 



P(x) = + 1 = {x + l){x^ - a; + 1) 

We claim that for a suitable starting value S{1) such as 5(1) = 2 we have: 

• S{n) = for infinitely many values of n. 

• S{n) = ^ (n + 2) and (n^ + n + l) are prime numbers. 
Here a table with the first values 



Values of n such that S{n) =0 


n + 2 


+ n + 1 


6{n + 2)d{n'' + n + l) 


3 


5 


13 


1 


69 


71 


4831 


1 


299391 


299393 


89635270273 


1 



P{x) = {2x-3){x'^ -x + 1) 

We claim that for a suitable starting value 5(1) we have: 

• 5(n) = for infinitely many values of n. 

• 5(n) = ^ (2n - 1) and (n^ + n + 1) are prime numbers. 
Here a table with the first values for 5(1) = 1, 2, 5 
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S{1) 


n such that S{n) = 


2n- 1 


n'-^ + n + 1 


1 


2 


3 


7 




24 


47 


601 




24186 


48371 


584986783 


2 


3 


5 


13 




69 


137 


4831 




658657 


1317713 


434093205307 


5 


6 


11 


43 




414 


827 


171811 




141629682 


283259363 


20058966965050807 



And all the values for 2n — 1 and + n + 1 in the table are primes. 

Above exemples allow us to unify all our previous observations (except for 
the conjecture 5 which is a special case needing more thought to be generalized). 
This is the next conjecture. 

5.1 Conjecture 6 

Clearly something very general is working and one can imagine to state a deep 
conjecture. Suppose P{x) = Hjli Qji^) where P is polynomial with integer co- 
efficients (not all zeros). Suppose each Qj is irreducible and Qi{k), Q2{k), . . . , Qm{k) 
can simultanously be primes for large k. Then we claim there are infinitely many 
values of S{\) such that: 

• S{n) = for infinitely many values of n. 

• S{n) = => Qj{n + 1) is simultaneously prime for 1 < j < m. 

5.1.1 RemEirk 

This could have direct application such as generating not only big primes (see 
conjecture 4) but also big primes with given property like twin primes. For 

instance let: 

• P{x + 1) = {x^ + l){x^ + 3) 
. 5(1) =4 

Then we compute big twin primes compared to the number of iterations: 



Values of n such that S{n) =0 


+ 1 


+3 


(5(n^ + l)S{n^ + 3) 


2 


5 


7 


1 


14 


197 


199 


1 


32374 


1048075877 


1048075879 


1 
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5.2 Conjecture 7 

Instead of considering an infinite sequence we consider a starting value and see 
what happens. This gives a modified version of the conjecture 6 with more 
details suggesting some analytic study. 

Suppose that P{x) is a polynomial with integer coefficients satisfying P{n) > 
and P(+oo) = +00. Let us define the sequence a as follows. 

• a(l) = N €N 

• a{n) = a{n — 1) — gcd(a(n — 1), P{n)) 

Then we conjecture without any other assumption on P: 

• 3 f{N) €{1,2,...,N} such that a{f{N)) = 0. 

• liminfjv^.(x) > and limsup n^oo^^ = 1- 

Now suppose as above P{x) = YVjLi Qji^) where P is polynomial with integer 
coefficients (not all zeros). Suppose each Qj is irreducible and Qi{k). Q2{k). . . . .Q. 
can be simultanously prime for large k. Then there are infinitely many N such 
that we get: 

• Qj{f{N) + 1) is prime for j G {1, 2, . . . , m} 

And more precisely there is a positive proportion of N such that Qj{f{N) + 1) 
is prime for j e {1,2,..., m} . Let us define this proportion as follows: 

. L(P) =lim„^^i#{fc I l<fc<n&nr=i'5(Q,(/(fc) + l)) = l}- 
Then we claim: 

• Qj is of degree 1 for any j G {1,2, . . . , m} => L{P) = 1. 

• One Qj is irreducible of degree 2 L{P) < 1 (see 3.3. with a sample 

giving L(2a;2 - 1) - 0.8). 

Here the graph of ^d{f{N) + 1) when P{x) = x foi 1 < N < 20000 

(fig-7) 
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Here the graph of i^S{2f{N) + 3)6{f{N) + 3)S{f{N) + 7) when P{x) = 
(2x + l){x + 2)(x + 6) for N = lOOfc and 1 < k < 20000 

(fig-8) 
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One can see the values of TV such that S{2f{N)+3)S{f{N)+3)6{f{aN)+7) = 
become sparse when N increases. We beheve that for N large enough this set 
is of measure zero since Qj are of degree 1. 



6 Conjecture 8 

We propose a way to prove again there are infinitely many twin primes but in a 
slightly different way than before. The recursion is the backwards version of the 
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Shevelev recursion. First we give a conjectural way to prove there are infinitely 
many prime numbers (this is very similar to what we said in conjecture 5 and 
should be easy to prove). Although it is not the simplest way to prove there 
are infinitely many primes, it seems important to start with the usual primes 
in order to see better how it is also working for twin primes. The two cases are 
indeed very similar. 

6.1 There are infinitely many primes 

Let: 

• a(l) = N-2>0 

• n>2 =>a{n) = a{n — 1) — gcd(a(n — 1), n — 1) 
Then we claim we have: 

1. \/N > 4, 3f{N) e{2,...,N-l} such that a{f{N)) = 0. 

2. f{N) ^ N {N ^ oo) and more precisely we claim /(TV) = N+oiN^^"^ log N) 

3. f{N) is prime. 

4. f{N) = N — 1 N — 1 is a. prime number. 

5. f{N) = N — 2 'i^ N — 2 is a,n odd prime number. 

6. f{N) = N — 3oiN — 4^N — 3 (resp iV — 4) is a prime number of form 
6k + 1. 

7. f{N) = N — 5oiN— 6^N — 5 (resp — 6) is a prime number of form 
30fc + 1. 

Since f{N) — >■ cxd as TV — >• c» (from 2.) there axe infinitely many primes (from 
3.). 

Here the graph of ^6{f{N)) ioi 2 < N < 20000. 

(fig-9) 
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To see the more precise behaviour claimed in 2. we plot — [jr^ for 3 < 
N < 20000. 

(fig.lO) 
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This is quite erratic and divided by log n it should converge very slowly to 
zero but more experiments are needed to confirm the tendancy. 

Remark 

It is worth to compare this sequence f{N) to p{tt{N)) where p{n) denotes the 
n-th prime and 7r(x) is the prime couting function. So that p{7t{N)) is the 
largest prime < N. Here we still plot ^-p(j(^)) for 3 < iV < 20000. 
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(fig- 11) 
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Things are less erratic and here the graph goes more certainly to zero without 
dividing by log n. We will discuss about an important consequence of this 
conjectured behaviour in section 9. 



6.2 There are infinitely many twin primes 

Let: 

• a(l) = A^ - 2 > 

• a(n) — a{n — 1) — gcd(a(n — 1), n + (—1)") 
Then we claim we have: 

1. V7V > 2, 3h{N) e {1, 2, . . . , A^} such that a{h{N)) = 0. 

2. h{N) - N {N ~^oo). 

3. VA^ > 98, h{N) and h{N) + 2 are necessarily simultanously primes. 

4. For iV > 4 we have h{N) = iV - 1 ^ (TV - 1, TV + 1) is a pair of twin 
primes. 

5. For iV > 13 we have h{N) — N — 2 4^ N is the greater of a pair of twin 
primes. 

6. For > 14 we have h{N) = N~34^N — 3is the lesser of a pair of twin 
primes. 

Since h{N) — >■ oo as A^ — > cx) (from 2) there are infinitely many twin primes 
(from 3). 

Here the graph of '^S{h{N))6{h{N) + 2) for 1 < Af < 30000 
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(fig- 12) 
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We see there is no more zero value for N > 98. 

7 Conjecture 9 

Let m e N and define the sequence a as follows: 

• a(l) e N 

• a(n) = \a{n — 1) — gcd(a(n — 1), mn + b{n))\ 

Where (6„)„>i is a periodic sequence of period length (3 such that: 

• mn + 6i, mn + b2, ■ ■ ■ , mn + 6^ can be simultanously prime for n large. 
Then we claim we have: 

• a(n) — for infinitely many values of n. 

• for n large enough we have a{n) = m{n + 1) + bi,m(n + 1) + 
62, ... , m(n + 1) + 6^ are simultanously prime. 

• the rate of growth of the sequence of nj such that a{ni) = is like (m + 1)' 

See APPENDIX 6 for experiments where we provide examples for /3-periodic 
sequences and /3 G {2, 3, 4, 5, 6}. 
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Generalisation 



Although we check few cases of this generalisation it is worth to mention this 
result which is an interesting variation of the conjectures 6 and 9. Suppose Qj is 
irreducible and Qi{k), Q2{k)^ . . . , Qp{k) can be simultanously prime for large k. 
Suppose (6n)n>i is a periodic sequence of period length /? such that {fe(«)}i<i<^ 
is a permutation of {*}i<i<^ define the sequence a as follows: 

• a(l) e N 

• a{n) = |a(n- 1) - gcd(a(n - l),Qb(n){n))\ 
Then we claim we have: 

• a{n) = for infinitely many values of n. 

• for n large enough we have a(n) = (3i(n + l), Q2{n + 1), . . . , Qfi{n^l) 
are simultanously prime. 

Extension of the conjecture 6 

We can extend what we done in conjecture 6 (and also in conjecture 7) with 
twin primes to prime triplet or any sort of m-uplet. For instance let us see how 
this is also working with prime triplet of type (p,p + 2,p + 6). Let (fe(f^))„>x be 
the 3-periodic sequence {2,6,0} and define the sequence: 

• a{n) = iV - 6 > 

• a(n) = a{n — 1) — gcd(a(n — 1), n + b(n)) 
Then we claim: 

1. V7V > 6 ^f{N) e {l,2,...,iV} such that a{f{N)) = 0. 

2. f{N) - N {N ^oo). 

3. V7V > 2735, (/(TV) + 1, f{N) + 3, f{N) + 7) is a prime triplet. 

4. For iV > 5 we have h{N) = N - I {N,N + 2,N + &) \s a. prime triplet 
of type {p,p + 2,p + 6) (sequence A022004 in ) 

Here the graph of ^S{f{N) + l)S{f{N) + 3)6{f{N) + 7) for 6 < Af < 30000 

(fig. 13) 
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We see that for N > 2735 there is no more zero value. 



8 Conjecture 10: on the Goldbach conjecture 

Finally we found a way to adapt the method to the Goldbach conjecture. This 
was not easy since this recursion is very sensitive (to the initial value and to 
what we put in the gcd). This formulation of the conjecture is nicer than the 
conjecture 5 and is very similar to the conjecture 6 and our formulation of the 
weak twin prime conjecture. Indeed the method works in both cases for N large 
enough. Namely define the sequence a for > 2 by: 

• ai = — 2 and for n > 2 by a„ — a„_i — gcd(a„_i, iV — {—1)^{N — n)) 
Then we claim: 

• there is a least G {2, 3, — 1} such that ag„ = . 

• grf ^ N {N ^ oo) 

• for iV > 2208 we have + 1 and 2N — g^ ~ 1 which are simultanously 
primes. 

Thus the Goldbach conjecture would be true. 

Hereafter a graph of ^S{gN + 1)S{2N - gN - 1) for 2 < TV < 30000. 

(fig. 14) 
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We see there is no more zero value for > 2208. 

Next we provide a zoom of this picture in the range 100000 < N < 130000. 

(fig. 15) 
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The Hues appearing regularly are due to the following observation. 

Relation with the weak twin prime conjecture 

Moreover we claim that Vm > 2 G N 

• dN — N — m for infinitely many values of A^. 
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And we add withouth any other condition on N: 

• gN ^N — m4^N — m+l and N + m — 1 are both primes. 
Thus for m = 2 this means there are infinitely many twin primes. 

The Goldbach constant 

It is worth to consider analytic aspects of this formulation of the Goldbach 
conjecture. For instance we claim: 

• ELs (l - t^) - C!y^i {n oo) where C < 4. 
Hereafter the graph of ti^i/^ Y,"^^^ (^1 - for l<n< 2000000 

(fig-16) 
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9 On a conjecture of Legendre 

In this study we give new formulation using gcd-algorithms for 3 problems 
among a list of 4 problems considered by Landau in 1912 as "unattackable 
at the present state of science" |Lan| |Pin| : 

• are they infinitely many primes numbers of form + 1? 

• the twin prime conjecture. 

• the Goldbach conjecture. 

So a Landau problem is missing but our gcd-algorithm formulation still works 
for this fourth probem attributed to Legendre. This remaining conjecture says: 
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• there is always a prime between 2 consecutive squares. 

Our belief is that for N large enough the algorithm given in 6.1 starting with 
a(l) = (TV + 1)^ — 2 produces a prime value f{{N + 1)^) greater than N'^. But 
this comes not from our asymptotic claim in 6.1. 

Remark 

Of course the situtation is clearer if we use the sequence p{7r{N)) introduced 
in 6.1. and the conjectured asmptotic behaviour which works for proving the 
Legendre conjecture (the behaviour toward zero of ^~p(^^)^ shown in 6.1. ). 
But for this one one must consider RH or things like that. Our hope is that the 
sequence f{N) despite its random nature allow us to avoid such considerations. 

A Goldbach variation 

We came also across a stronger hypothesis than the original Legendre conjecture 
and it is worth to present this "Goldbach variation". We don't know if this last 
conjecture could be helpful for finding a "gcd-formulation" of Legendre conjec- 
ture like our conjecture 10 for the Goldbach conjecture (we found nothing clear 
even for N large enough). Namely we claim: 

• V7V > 175 there is at least one k satisfying 2 <k <2N such that N'^ + k+1 
and {N + 1)^ — k are simultanously primes. 

Here the graph of the function of for 2 < iV < 20000 representing the number 
of k verifying 2 < k < 2N and such that N'^ + k + 1 and {N + 1)^ - fc are 
simultanously primes. 

(fig- 17) 



450 
A-00 




It looks like the classical Goldbach comet. 
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Concluding remark 

The Schinzel hypothesis doesn't cover the Goldbach conjecture but in view of 
our conjecture 10 we think there is a possible extension of the conjecture 6. We 
mention that some people didn't think Rowland simple theorem and related 
ideas of V. Shevelev could lead to interesting results regarding the theory of 
the prime numbers. We hope we show this advice is wrong and that some very 
important rules need to be better understood and that our results could inspired 
some analytic or probabilistic studies. Regarding the probabilistic approach 
we think we have other striking arguments. In our study we consider only 
natural objects in the gcd such as polynomials or periodic sequences. It appears 
we can play with less predictive functions and for instance we observe we can 
generate twin primes with any function r taking values in {0, 2} sufficiently 
"randomly"!^ To see this let us consider the differences of the Beatty sequence 
for TT i.e. Vn = \T^n\ — \n{n— 1)J which takes values in {3,4} and is not a 
periodic sequence (this is the sequence A063438 in [Slo' ) . 

Then let r„ — 2(i'„ — 3) wich takes values in {0, 2} and define the sequence: 

• a(l) = TV - 2 > 

• a{n) — a{n — 1) — gcd(a(n — 1), n + r„) 
Then we claim: 

1. V7V > 2, 3/(7V) € {1, 2, . . . , iV} such that a{f{N)) = 0. 

2. f{N) - N {N ^oo). 

3. VA^ > 1649, f{N) + 1 and f(N) + 3 are necessarily simultanously primes. 
Hereafter we plot ^(5(/(A) + l)(5(/(iV) + 3) for 2 < iV < 30000. 

(fig.18) 

^This needs to be clarified but we think lim„_>cx> ^ 5Zj=i ^ {^i 2} is a sufficient condition 
and we believe the random process comes mainly from the gcd. 
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We can see there is no more zero value for > 1649. 
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APPENDIX 1 



Table related to the conjecture 1 



ft- 








1 


2 


1 


1.000000000 


2 


5 


1 


1.250000000 


3 


11 


1 


1.375000000 


4 


23 


1 


1.437-")00000 


5 


47 


1 


1.468750000 


6 


79 


1 


1.234375000 


7 


157 


1 


1.226562500 


8 


313 


1 


1.222656250 


9 


619 


1 


1.208984375 


10 


1237 


1 


1.208007812 


11 


2473 


1 


1.207519531 


12 


4909 


1 


1.198486328 


13 


9817 


1 


1.198364257 


14 


19603 


1 


1.196472167 


15 


39199 


1 


1.196258544 


16 


''70 1 no 

78193 


1 


1.193130493 


17 


156019 


1 


1.190330505 


18 


311347 


1 


1.187694549 


19 


622669 


1 


1.187646865 


20 


124414!) 


1 


1.18(r)12947 


21 


2487739 


1 


1.186246395 


22 


4975111 


1 


1.186158895 


23 


9950221 


1 


1.186158776 


24 


19900399 


1 


1.186156213 


25 


39800797 


1 


1.186156183 


26 


79601461 


1 


1.186154201 


27 


159202369 


1 


1.186150081 


28 


318404629 


1 


1.186149675 


29 


636788881 


1 


1.186111720 
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APPENDIX 2 



In the following table one can see the number of iterations n giving a{n) = 
(loft column) and the corresponding value mn + m — 1 in the middle column. 
The right colum gives the value of the function ''isprime" using pari-gp. 



TO = 5 



n such that a{n) = 


5n + 4 


S{5n + 4) 


2 


14 





17 


89 


1 


95 


479 


1 


575 


2879 


1 


3419 


17099 


1 


19967 


99839 


1 


119801 


599009 


1 


718571 


3592859 


1 


4311419 


21557099 


1 


25867229 


129336149 


1 


TO = 6 


n such that a(n) = 


6n + 5 


6{6n + 5) 


2 


17 


1 


16 


101 


1 


7() 


4()1 


1 


466 


2801 


1 


3258 


19553 


1 


22774 


136649 


1 


159306 


955841 


1 


1114124 


6684749 


1 


77796204 


46777229 


1 


54573434 


327440609 


1 


TO = 7 


n such that a(n) = 


7n + 6 


5{7n + 6) 


2 


21 





23 


167 


1 


113 


797 


1 


899 


6299 


1 


6973 


48817 


1 


55633 


389437 


1 


444901 


3114313 


1 


3558575 


24910031 


1 


28468585 


199280101 


1 
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m = 8 



n such that o(n) = 


8n + 7 


S{8n + 7) 


2 


23 


1 


20 


167 


1 


188 


1511 


1 


1682 


13463 


1 


15020 


120167 


1 


134504 


1076039 


1 


1210544 


9684359 


1 


10894874 


87158999 


1 


98053784 


784430279 


1 


m = 9 


n such that a{n) = 


9n + 8 


6{9n + 8) 


2 


26 





29 


269 


1 


299 


2699 


1 


2935 


26423 


1 


28869 


259829 


1 


288;-')8.") 


2.19.147;-! 


1 


2883809 


25954289 


1 


28832339 


259491059 


1 



m = 100 



n such that a{n) — 


lOOn + 99 


d{100n + 99) 


2 


299 





226 


22699 


1 


22810 


2281099 


1 


2303908 


230390899 


1 



To see further how this is working nicely let us take m = 1000. 



n such that a(n) = 


lOOOn + 999 


^(lOOOn + 999) 


2 


2999 


1 


2986 


2986999 


1 


2917174 


2917174999 


1 



Finding big primes (from 1.3.1) Wc take m = 2^ for 100 < k < 200 so 
that b{l) = 3 and a{b{l) + 1) = 3.2'^ — 1 is the first record value (not necessarily 
prime of course and it would be interesting to find conditions on the initial value 
in order to have very few terms to compute) . Next we stop the algorithm after ~ 
10^ iterations (few seconds are needed each time). In the following table we keep 
the values of k such that {rn+l)a{b{l) + l)+m+mJ2j:^4 {<^{j + 1) ~ f^O) + 1) 
is a prime value of size ~ 3.4'^and should be our second record value a{b{2) + 1). 
We notice we came across 11 primes in that range. 
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k 


a(&(2) + 1) 


100 


4820814132776970826625886270541990288599672051495735016816639 


107 


78984218751417890023438520762741628349070105828337814558023352319 


118 


331283824645947061796868281389238893531663193824530694464237235871940607 


127 


86844066927987146567678238756515930889442064948849015334398126094787194912767 


131 


22232081133564709521325629121668078276785918415807297904592970941163697569005567 


132 


88928324534258838085302516486672313230587227346404196132053363353462480809492479 


141 


23312026706708748851033542881882226879708773352492831418233933930240272261700160323583 


149 


1527776982250864564701334266307033620788600839243814368422923338893418735708806455350001663 


158 


400497569235170640449066569906791021488007090321237481912409760223102559592991367888359129087999 


164 


1640438043587258943279376670338216024014877072467663636588378426222267002018688816605037919597494271 


172 


107507747624534602106757229467285325349838983867263171688888770338883656854354650572386871166337445527551 



pari gp code 

f or (k=100 , 200 ,m=2-k ; a=l ; S=0 ; M=0 ; 

for(n=2,10-5,t=a;a=abs(a-gcd(a,m*n-l)) ;M=M+lf (a-t>0,a-t,0) ;S=S+lf (a-t<0,a-t+l,0) ; 
if (abs(t-a)>l,if (isprlme((m+l)*M+m+m*S)==l,priiit(k," ", (m+l)*M+m+m*S, "")))) ) 
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APPENDIX 3 



Let a(l) = m and 

• a{n) = \a(n — 1) — gcd(a(n — 1), — 1)| . 

Here a table for various starting values m allowing us to exhibit distinct pairs 
of twin primes (the other values m = lOA; for 1 < A; < 14 produce also pairs of 
twin primes but they are all in this list keeping the distinct pairs only). 



Tfl 


Ti snrVi tli/^t nln) = 

IV O lAvyXJ. UXJ-OiU LV\lljf \J 


Sin) 


Sin + 2") 




1 1 
ii 


i 


1 

i 




1 0'7 

io7 


i 


1 






1 

i 


1 


zL) 


1 T 

17 


i 


-l 

1 






i 


1 

i 




79009 


i 


1 

i 




29 


i 


1 

1 






i 


1 

1 




777Uii 


i 


1 

i 


/in 


/1 1 


1 
i 


1 
i 




1787 


1 


1 




3198731 


1 


1 


60 


59 


1 


1 




3527 


1 


1 




12448001 


1 


1 


70 


71 


1 


1 




5099 


1 


1 




26010041 


1 


1 


100 


101 


1 


1 




10499 


1 


1 




110258891 


1 


1 


110 


107 


1 


1 




11699 


1 


1 




136890881 


1 


1 


140 


137 


1 


1 




19181 


1 


1 




367953497 


1 


1 



And 2 more starting values 



m 


n such that a{n) = 


S{n) 


S{n + 2) 


200 


197 


1 


1 




39161 


1 


1 




1533646397 


1 


1 


300 


281 


1 


1 




79559 


1 


1 




6329815697 


1 


1 



36 



We also provide this computation but using Rowland-Shevelev recursion: 

• ai = 3 then a„ = a„_i + gcd(a„_i,n(n — 2)) 

Here the table of records of the differences a„ — an-i which are the lower terms 
of some twin prime pairs. 

records of a„ — a„_i 
3 

5 

n 

41 ~ 

101 

239 

521 

1049 

2111 

4229 ~ 
10331 

20747 ~ 

41519 

83219 

166847 

333791 

669479 

1341017 ~ 

2682539 

5365229 

10732751 

21466259 

42932567 

85865321 

171730679 

343461647 

686929511 

1373891861 

2747784329 ~ 

5495586839 
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APPENDIX 4 



Here a(l) = 4m^ and a{n) = \a{n — 1) — gcd(a(n — 1), n(n + 2m))| and the 
3 first values of n such that a„ = . 



lit 


— 1— 1 TKTr\ ATI ni Ti\ — n 

/fc 1^ J- WllCll \Ai\lb\ U 




R{ri A- -\- ^^ 







1 
i 


1 
1 




193 


1 


1 




38197 


1 


1 


3 


31 


1 


1 




1091 


1 


1 




1197193 


1 


1 


4 


59 


1 


1 




4013 


1 


1 




16138511 


1 


1 


5 


97 


1 


1 




10477 


1 


1 




109880317 


1 


1 


6 


139 


1 


1 




20521 


1 


1 




421370778 


1 


1 
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APPENDIX 5 



p = 2 



Values of n such that a(n) = 


2(n + 1)2 - 1 


5(2(n+ 1)2-1) 


3 


31 


1 


35 


2591 


1 


2627 


13812767 


1 




287()8012()871111 


1 


p = 3 


Values of n such that a{n) = 


3(n+ 1)2 - 1 


(5(3(71+1)2-1) 


3 


47 


1 


51 


8111 


1 


7665 


176302667 


1 


176310323 


93255991046954927 


1 



p = 5 



Values of n such that a{n) =0 


5(n+ 1)2 - 1 


(5(5(n+ 1)2 - 1) 


3 


79 


1 


83 


35279 


1 


34647 


6002419519 


1 



p = 7 



Values of n such that a(n) =0 


7(n+ 1)2 - 1 


<5(7(n+ 1)2-1) 


5 


251 


1 


257 


465947 


1 


461009 


1487711540699 


1 



p = n 



Values of n such that a(n) = 


ll(n + l)2 - 1 


^(ll(n + l)2 - 1) 


11 


1583 


1 


1419 


22180399 


1 


22181509 


5412213244681099 


1 


p=l7 


Values of n such that a{n) = 


17(n + 1)2 - 1 


(5(17(n + 1)2 - 1) 


11 


2447 


1 


2417 


99394307 


1 


87543523 


130285766103755791 


1 
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APPENDIX 6 

We define the sequence a as follows: 



• a{n) = \a{n — 1) — gcd(a(n — 1), mn + b{n))\ 

We then give tables supporting the conjecture 9 for various k, m and b periodic 
sequence which is given by its period {bi, 62, b^}. 

2 periodic sequence 

b = {0.2} . m = l.k = 100 



n such that a{n) — 


d{n+ l)(5(n + 3) 


100 


1 


196 


1 


310 


1 


616 


1 


1228 


1 


2380 


1 


4648 


1 


8860 


1 


17()2(i 


1 


33808 


1 


67408 


1 


134680 


1 


267718 


1 


535348 


1 


1069216 


1 


2138398 


1 


4275640 


1 


8545696 


1 


17091376 


1 


34182748 


1 


68365468 


1 


136730638 


1 


273461158 


1 


546917140 


1 


1093813726 


1 


2187610990 


1 



40 



6={-l,l},m = 10,fc = 100 



n such that a(n) =0 


<5(10n + 9)^(10n + 11) 


101 


1 


1115 


1 


12203 


1 


130013 


1 


1427183 


1 


15692309 


1 


172614683 


1 



b = {-2,2} ,m = 3,k = 100 



n such that a{n) = 


5(3n + l)(5(3n + 5) 


34 


1 


116 


1 


434 


1 


1576 


1 


6102 


1 


21154 


1 


84606 


1 


338386 


1 


1351382 


1 


5405526 


1 


21622094 


1 



3 periodic sequence 

6 = {2,6,0},m = l,fc = 3000 
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n such that o(n) =0 


d{n + l)6{n + 3)6{n + 7) 


2686 


1 


5230 


1 


10456 


1 


19420 


1 


29566 


1 


54496 


1 


105526 


1 


211060 


1 


408430 


1 


802126 


1 


1600216 


1 


;-')20()2()() 


1 


6393910 


1 


12783496 


1 


25566676 


1 


51095410 


1 


102190390 


1 


204347176 


1 


6 = {0,4,6},TO = l,fc = 3000 


n such that a{n) = 


S{n+l)6{n + 5)6{n + 7) 


2082 


1 


3462 


1 


6546 


1 


12372 


1 


23052 


1 


44262 


1 


85086 


1 


167016 


1 


313986 


1 


622476 


1 


1237206 


1 


2452752 


1 


4882326 


1 


9753276 


1 


19504866 


1 



b = {-2, 2, -4} , m = 5, A: = 2000 
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n such that a{n) = 


6{5n + l)(5(5n + 3)S{5n + 7) 


1772 


1 


9806 


1 


58274 


1 


343772 


1 


2057378 


1 


12342518 


1 


73895180 


1 
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4 periodic sequence 

b = {1, 7, 11, -1} , m = 1, A; = 20000 



n snrb that n,(7i?\ = 


S(n)S(n + 2)S(n + 8)6(n + 12) 


1 O/l 01 


1 

i 


6o( (\) 


1 

i 


/Uo4i 


-1 
i 


138239 


1 


236771 


1 


443159 


1 


882239 


1 


1758389 


1 


3376979 


1 


6631901 


1 


13236539 


1 


26425379 


1 


52658999 


1 


104785649 


1 


209560319 


1 


418973999 


1 



b = {2, 6, 18, 26} , m = 1, fc = 100000 



n such that a{n) — 


^•=l'^(r^ + Ki) + l) 


83200 


1 


150190 


1 


294754 


1 


573844 


1 


1107784 


1 


22()8()()4 


1 


4171774 


1 


8332840 


1 


16461094 


1 


32756680 


1 


65166814 


1 


130175344 


1 


260331034 


1 


520484380 


1 


1040389234 


1 


2080515244 


1 


4161006904 


1 


8321226490 


1 
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5 periodic sequence 

b = {2, 8, 12, 14, 18} , m = 1, fc = 1500000 



n such that a{n) — 


TT'' 1 d(n + l + b(j)) 

LL'j — l \ ' ' \J / / 


1212424 


1 


2270674 


1 


4271158 


1 


8358658 


1 


15875398 


1 


31562608 


1 


62555878 


1 


125087098 


1 


249509788 


1 


477331018 


1 


954642034 


1 


1905881278 


1 


3809937208 


1 


h = {2. 8. 12. 14. 18} . rn = 1. k = 2000000 


n such that a{n) = 


n-=i'5(n + l + Hi)) 


1460728 


1 


2839924 


1 


4218154 


1 


8068438 


1 


16130884 


1 


32240278 


1 


64123234 


1 


127725328 


1 


254416288 


1 


507764278 


1 



6 periodic sequence 

b = {0, 2, 6, 14, 30, 62} , m = 1, fc = 2000000 
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n such that a{n) = 


n-=i'^(r^ + l + &(i)) 


1460728 


1 


2839924 


1 


4218154 


1 


8068438 


1 


16130884 


1 


32240278 


1 


64123234 


1 


127725328 


1 


254416288 


1 


507764278 


1 
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